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CN ■ Abstract 
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' The polarization of radiation by scattering on an atom embedded in combined exter- 

Qh| nal quadrupole electric and uniform magnetic fields is studied theoretically. Limiting 

O 

' cases of scattering under Zeeman effect and Hanle effect in weak magnetic fields are 

■ discussed. The theory is general enough to handle scattering in intermediate mag- 

■ netic fields (Hanle-Zeeman effect) and for arbitrary orientation of magnetic field. 
^ . The quadrupolar electric field produces asymmetric line shifts and causes interesting 

level-crossing phenomena either in the absence of an ambient magnetic field or in 
its presence. It is shown that the quadrupolar electric field produces an additional 
depolarization in the Q/I profiles and rotation of the plane of polarization in the 
U 1 1 profile over and above that arising from magnetic field itself. This characteris- 
tic may have a diagnostic potential to detect steady state and time varying electric 
fields that surround radiating atoms in Solar atmospheric layers. 
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1 Introduction 



Scattering of polarized radiation by an atom is a topic of considerable interest 
to astrophysics ever since Hale [T] first observed polarization related to Zeeman 
effect in spectral lines originating in Sun spots. The polarized radiation is usu- 
ally expressed in terms of the Stokes parameters. The concept of scattering ma- 
trix connecting the Stokes vector S' of incident radiation to the Stokes vector 
S of scattered radiation was introduced quite early in the context of Rayleigh 
scattering [2] . Polarized radiation in spectral lines formed in the presence of an 
external magnetic field has been studied widely and a comprehensive theoreti- 
cal framework has been developed [3pil5]|6ll71l8]|9lll0llllfl2fl3lll4fl5j . The Hanle 
effect is a depolarizing phenomenon which arises due to 'partially overlapping' 
magnetic substates in the presence of weak magnetic fields, when the splitting 
produced is of the same order or less than the natural widths. Favati et al. [16] 
proposed the name 'second Hanle effect' for a similar effect in 'electrostatic 
fields'. Casini and Landi Degl'Innocenti [17j have discussed the problem in 
the presence of electric and magnetic fields for the particular case of hydrogen 
Lyman a line. It was followed by a more recent paper by Casini [18j. The 
relative contributions of static external electric fields, motional electric fields 
and magnetic fields in the case of hydrogen Balmer lines, have been studied by 
Brillant et al. [19]. A historical perspective and extensive references to earlier 
literature on polarized line scattering can be found in Stenfio [12j, Trujillo 
Bueno et al. [20] and Landi Degl'Innocenti and Landolfi |21j . 

A quantum electrodynamic theory of Hanle-Zeeman redistribution matrices 
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has been developed by Bommier [TOpT] and Landi Degl'Innocenti &; Co- 
workers (see the book by [21j). The formulation presented in [TO] and [H] 
includes the effects of partial frequency redistribution (PRD) in line scatter- 
ing for a two-level atom. It is a perturbation theory, in which PRD effects 
appear in the fourth order. The theory presented in [21] and references cited 
therein, considers only complete frequency redistribution (CRD) in line scat- 
tering. 

A classical theory of line scattering PRD for the Hanle-Zeeman effect has been 
formulated by Bommier & Stenfio [15j. This theory is non-perturbative and 
describes the scattering process in a transparent way. The classical theory for 
Hanle-Zeeman scattering developed by Stenfio [2] considered only coherent 
scattering in the laboratory frame. In [TJ] the redistribution matrices were 
derived in the atomic rest frame. The corresponding laboratory frame redis- 
tribution matrices have been derived in [22] • The equivalence between the 
classical (non-perturbative theory) and quantum electrodynamic (perturba- 
tive theory) redistribution matrices for the triplet case is established in [231. 
In all these papers only the dipole type line scattering transitions in the pres- 
ence of pure magnetic fields is considered. Taking into account all higher order 
multipoles as well, polarization of line radiation in the presence of external 
electric quadrupole and uniform magnetic fields was studied [241125] . where 
scattering of radiation by atoms, however, was not considered. 

The purpose of the present paper is to develop a quantum electrodynamical 
approach to scattering processes in the presence of external electric and mag- 
netic fields of 'arbitrary strengths', taking also into consideration all other 
multipole type transitions apart from the usually dominant electric dipole 
transition. The atomic electron is represented using non-relativistic quantum 
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theory including spin. The radiation field is described in terms of its electric 
and magnetic multipole states, in a second quantized formalism. The external 
electric field is assumed to be 'quadrupolar' in nature, while the magnetic field 
is uniform and arbitrarily oriented with reference to the principal axes frame 
of the electric quadrupole field. This general formalism can be employed also 
to solve the scattering problems involving linear steady state electric fields at 
the radiating atom. 

In Sect. 2 we describe the theoretical formulation. In Sect. 3 the scattering 
matrix for the general physical situation is derived. The particular case of the 
dipole transitions for a triplet is also considered, for the purpose of comparison 
with Stenfio [H] in the pure magnetic field limit. Sect. 4 contains numerical 
results and discussions. Conclusions are presented in Sect. 5. 

2 Theoretical formalism for scattering 

We consider polarized radiation incident on an atom along an arbitrary di- 
rection {6' , (j)') and getting scattered into a direction {6, (p) with respect to a 
conveniently chosen right handed Cartesian coordinate system, referred to as 
the Astrophysical Reference Frame (ARF) and shown as (X, F, Z) in Fig. [H If 
v' and V denote respectively the frequencies of the incident and scattered radi- 
ation, we may define wave vectors k' and k with polar co-ordinates {k\ 6', (f)') 
and (fc, 9, (p) where k' = 2'kv' = uj' and k = 2ttv = uj in natural units with 
h= 1, c = 1 and mass of the electron me = 1. The atom is exposed to an exter- 
nal magnetic field B with strength B, directed along {6b,4>b) and an electric 
quadrupole field characterized by strength A and asymmetry parameter t] in 
its Principal Axes Frame (PAF), which is denoted by {Xq, Yq, Zq) in Fig. [H 
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The transformation to PAF from ARE is achieved by a rotation R{aQ, /3q, 7q) 
through Euler angles (aQ,/5Q,7Q) as defined by Rose [26]. The magnetic field 
B is directed along {6b, (Pb) with respect to PAE. Eollowing Rose [26j, we 
define left and right circular states of polarization {fi = ±1) respectively, 
which are mutually orthogonal to each other and to k. We use here the symbol 
s^=±i instead of Up=±i employed in Rose PB]. Like wise e'^,^_^^, which are or- 
thogonal to k'. Any arbitrary state of polarization e of the incident radiation 
may then be expressed as e = c_^_i e_^_i + c_^ e_i using appropriate coefficients 
c'_j_^, which are in general complex and satisfy |c'|_]^p + |c'_^p = 1. We, therefore, 
denote the orthonormal states of polarized incident radiation by |k',/i'), with 
/i' = ±1. We seek the probability for scattering into two polarized states of 
scattered radiation |k, /i), /i = ±1 on an atom which is initially in a state ipi 
with energy Ei before scattering and makes a transition to a final state ipf 
with energy Ef, in the process of scattering of polarized radiation. 

2.1 Energy levels of an atom in electric quadrupole and uniform magnetic 
fields 

The energy levels of an electron in an atom are primarily determined by the 
Hamiltonian 

Ho = -\ +V{r) , (1) 
where V{r) denotes its Coulomb interaction with the nucleus. If we start with 
the Dirac equation [17J and use its non-relativistic reduction, terms like spin- 
orbit interaction may also be included in Hq. In the absence of external fields, 
the energy levels of the atom are determined by 

H^ = ^:Hoi^)+ j: (2) 

i=l i>j=l *i 
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where e denotes the charge of the electron, = |rj — rj| and Z denotes the 
atomic number. If E denotes an energy level and tl) the corresponding wave 
function of the atom with total angular momentum J, it is well known in 
the context of Zeeman effect that E gets split into (2 J + 1) equally spaced 
levels Em = E + gBM with corresponding energy eigenstates | JM)b, M = 
J, J — 1, ■ ■ ■ , — J + 1, — J, when the atom is exposed to an external uniform 
magnetic field B with strength B. The states |JM)b are defined with the 
axis of quantization chosen along B and g denotes the magnetic-gyro ratio or 
Lande (yf-factor. For B < 100 gauss, when gB is of the same order as the width 
of a line, Hanle effect [27j takes place. For a line in the optical range, the region 
100 < B < 1000 gauss is generally referred to as the Hanle-Zeeman regime. 
For B < 1 gauss, one has to pay attention to the interaction of electron with 
the magnetic and electric moments of the nucleus, which give rise to hyperfine 
splitting [281129] . If the atom is exposed to an external electric quadrupole field 
either by itself or in combination with B, the splitting of the energy levels is 
not, in general, equally spaced [30|31II24|[25] and in such scenarios, the atomic 
Hamiltonian in PAF is given by 



Ha = H^ + g 3 -B + A 



(3) 



The split energy levels may be denoted by Eg., where s takes values s = 
1, 2, ■ ■ ■ , (2J+ 1) starting from the lowest level (s = 1) for a given J. The cor- 
responding energy eigenstates may be denoted by \J,s), which are expressible 

as 

J 

\J,s)= J2 <M,V,B,9bAb)\JM)q, s = 1,2,---,(2J+1) (4) 

M=-J 

where \JM)q are defined with the quantization axis chosen along the Z-axis, 
Zq of the PAF. The notation c^^ was used in for J = 1,3/2 to denote 
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the expansion coefficients, without any specified convention for ordering of the 
levels. We may rewrite Eq. (jlj) as 

J 

= ^ 4 \ Jm) , (5) 

in terms of the \Jm) states, which are defined with respect to the Z-axis of 
ARF chosen as the quantization axis. Clearly, 

J 

M=-J 

and hence depend on ag, (3q, 7q, A, r], B. If the magnetic field is absent, the 
depend only on aq, (3q,'~^q, A,rj since a%j in that case [24j depend only on 
A and rj. It may be noted that the frame of reference employed in was 
PAF itself, i.e., aQj^Qj^jq = 0; hence 9b,4>b was used there instead of the 
9b, (pB here and the Euler angles {aq, Pqj'Jq) find no mention there. On the 
other hand, if the electric quadrupole field is absent and the atom is exposed 
only to a magnetic field B directed along {9b, (Pb), it is clear that 

C= E DU{<Pb,9b,0)6m,s-j-i, (7) 

M=-J 

which reduces to 

Cm = ^s,J+m+l , (8) 

if the field B is along the Z-axis of ARF itself. 

In general, therefore, when the energy levels of an atom are defined through 
Ha^Ju = Enipn, the atomic wave functions ipn are of the form given by Eq. ([5]), 
which specialize appropriately to \JM)b or \Jm) if Eq. ([7]) or Eq. are 
used instead of Eq. ([6]). Thus, in general, the complete set of orthonormal 
energy eigenstates of an atom in a combined external electric quadrupole and 
uniform magnetic field environment may be denoted by {ipn}, where n is used 
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as a collective index, which includes the serial number Sn along with the total 
angular momentum J„ and all other quantum numbers which may be needed 
to specify each ipn uniquely. In the presence of a pure magnetic field B, the 
magnetic quantum number Mn replaces s„ through the 5-function in Eq. ([7]). 
Moreover, if B is along the Z-axis of ARF itself, s„ gets replaced by m„ 
through the 5-function in Eq. ([H]). 

In general, a summation over n as in l'^n)(^n| = 1? implies a summation 
with respect to s„ as well. This summation over s„ may be replaced by a 
summation with respect to M„ or m„ in some particular CcLSGS cLS mentioned 
above. The initial and final states of the atom before and after scattering 
are denoted by ipi and ipf. They also belong to {ipn}- We use the short hand 
notation 

K) = |V'.;k',/i'); |/) = |^^;k,/x). (9) 

2.2 Interaction of atom with the radiation field 

It is well known that the local minimal coupling i.e, ip'^yipAy (with implied 
summation over v) of the Dirac field if) and the electromagnetic field rep- 
resented by the four potential A^, i/ = I,-- - ,4 is the fundamental interac- 
tion responsible for all electrodynamical process involving photons and elec- 
trons [32|33] . In the interaction representation, ip and Ay satisfy the free field 
equations of Dirac and Maxwell respectively. The quantity ip = ^^74, where 
ip"^ denotes the hermitian conjugate of ip and 71, 72, 73, 74 are 4x4 Dirac ma- 
trices. To facilitate calculations using the atomic wave functions {ipn}, we may 
use the non-relativistic two componental forms of ip and ip in c number theory 
for electrons, retain the Maxwell field in q number theory and represent the 
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interaction of the radiation field in tlie Coulomb gauge with the atom as 



-lAivj, t) • Vi + 2 ■ ( VjX A(rj, t) 



-iHAt 



(10) 



where aj denote the Pauli spin matrices of the electron labeled j located at 
Tj and Z denotes the atomic number. The quantum field variable A(r,t) in 
interaction representation may be expressed as 



+ A 
k"/i"^k"^t" 



at,,,,Airir^i/(r)*e^^ * 



'111 



where cu" = |k"| and the creation and annihilation operators, denoted by a^^ 
and ak/i respectively, satisfy the commutation relation 



5(k - k') 6^ 



'12) 



for any pair k, n and k', /i' in general, while 



Ak^(r) e^''-" , 



(13) 



denotes a c number and Ak^(r)* denotes its complex conjugate. In particular, 
the operators are also used to generate the initial and final states of radiation 
in Eq. ([9]) through 

|kV) = a+^,| )o; (k/i| = o( kk/. , (14) 
where | )o denotes the vacuum state of the radiation field. 



9 



2.3 The Scattering process 

The S'-matrix for scattering may be defined, as usual [34J, by 

5 = lim t^oo Uit.to) , (15) 

tQ — * — oo 

where the evolution operator satisfies 

U{t, to) = 1-1 f dt'Hint{t')U{t\ to) , (16) 
J to 

which on iteration leads to the perturbation series 

°° roo rti ftN-l 

^ = 1 + E i-^)"" / / d-^^--- dt^H,nt{ti) ■ ■ ■ Hi^tN) , (17) 

J~oo J—oo J —oo 

since Hint{t) given by Eq. flTU]) is linear in A (see Eq. f lTT]) ). the first order 
(A^ = 1) term can contribute to either absorption through the first term in 
Eq. (ITT!) or to emission through the second term in Eq. (ITT!) and the integral 
over dti, from — oo ^ oo leads to the respective energy conservation criteria 
of Bohr. In the scattering problem under consideration, the lowest order (in e) 
contribution to {f\S\i) is obtained from the N = 2 term, which we may denote 
as {f\S'^\i). We introduce Enl^n)(^n| = 1 between Hint{ti) and Hint{t2), 
neglect contribution from two photons in the intermediate state and employ 
the notation |n) = \'>pn)\ )o- This leads, on using Eqs. (fTTj) and (fT2|) . to 

{n\HUm) = A.(k',/i')e[^(^"-^-"')*^] , (18) 

{f\HUti)\n) = ^/„(k,/i)eW^/+--^")*i] , (19) 
where ^„j(k',/i') and £^/^„(k, /i) denote amplitudes for absorption and emis- 
sion, involving Ak'|j'(rj) and Ak^^^j)* respectively, which are independent of 
time variable, instead of A{r j,t). We may change the variable of integration 
from t2 to ^2 = t2 — ti, ranging from — oo 0, associate a width r„ with ip^ by 
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introducing a factor exp(Tnt2) (see [35)136] ) and obtain after completing both 
the integrations, the expression 

(/|^2|^) = -2 TT z 6{Ef + uj-Ei-uj') Tfi{ki^; k'/i') , (20) 

where the on-energy-shell T-matrix element is of the form 

Tfi(k,fi;k',fi') = ^^SfniKfi) (pn Xj(k',/i') , (21) 

n 

and the profile function is given by 



'^n= {^nf - i^n) ^ ^nf = En-Ef, (22) 



on making use of En — Ei — uj' = ujni — uj' = Unf — u hj virtue of the energy 
(5-function in Eq. (!20l) . Using Eq. (jlj) and observing that {ipn} are completely 
antisymmetric with respect to the labels 1, 2, ■ ■ ■ , j, ■ ■ ■ , Z of the electrons, we 
have 

m'„,m'- 

£/„(k,/.)= 4c„(eJX^|£:(k,/i)|J„m:), (23) 

In 

where the matrix elements on the right hand side satisfy 

(J„M„|^(k,/i)|J,mi) = {Jimi\£{k, ^i)\J^m^Y , (24) 

between any pair of lower and upper atomic states and 

Ze 



^(k,/i) 



(27r)3/2^ 



-^A^^^ . V + -a • (V X A*^^) 



(25) 



with respect to an electron in the atom. Since atomic transitions during ab- 
sorption and emission conserve total angular momentum and parity, we use 
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the standard multipole expansion [26| for Ak,;^ given by Eq. (fT3l) . viz, 

oo L 



^LAf I'll ±^ 



LM 



, (26) 



L=l Af =-L 

where [L] = {2L + 1)^/^ and (6*, 0) denote polar angles of k, while A^} and 
A^]^ denote respectively the 'magnetic' and 'electric' 2^-pole solutions of the 
free Maxwell equations. Using the notation 



27rV2cj 



-^A^^) ■ V + • (V X a£/^) 



(27) 



and noting that Eq. (1271) is an irreducible tensor of rank L, we may apply the 
Wigner-Eckart theorem to write 



{Jumu\A(k,fx)\Jimi) = A{k, (x)m^rni = YC{Ji, L, Ju;mi,M,mu)JLi 

L 

x(z^)^+(^)D^. (0,^,0), 



(28) 



where the reduced matrix elements are given by 



Moo) = {Ujt\uj)U)g4L) + {Ujt\u)U)gM , (29) 



in terms of the projection operators 



9±iL) 



1 ± (-1) Vvr, 



(30) 



In the above equation vr^, tt; denote the parities of the upper and lower levels. 
Using Eq. (^^, we have 



{Jimi\S(k,fi)\Jumu) = £(k,n)mimu=YCi'Ji,L, Ju.m, M,mu)jL\ 



00 



XI 



-^/x)^+(^)D^,,X0,^^,O)*. (31) 



Thus, we may express Eq. (pTi) as 



T/i(k;u;kV 



f ""i 



^(k,/i)e;^M(k',/i' 



(32) 
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where the summation J2n imphes summation with respect to s„ as well, and 
.4(k',/i') and £^(k, /i) denote matrices, whose elements 

(J„m'„ |^(k', fx')\Jim[) = Aik', fi')m'„m':, 

(J/m).|£:(k,/i)|j;„m^) =£:(k,/i)„'^^// , (33) 

may be written explicitly using Eqs. (128|) and (13T!) and Q denotes a hermitian 
{2Jn + 1) X {2Jn + 1) matrix, which is defined in terms of its elements 

Sm'^m'„ = Cm'J, • (34) 

Clearly, the summation over n on right hand side of Eq. (132|) indicates a sum- 
mation with respect to all the atomic states {ipn}, which constitute the com- 
plete orthogonal set. Since n is a cumulative index Y.n includes Y.j„ I]s^=l^) 
apart from summation with respect to other quantum numbers. The left hand 
side of Eq. fl32|) is written for a given ipi and tpf with energies Ei and Ef re- 
spectively. The quantities Si and s / are specified by left hand side of Eq. (1521) 
and hence they are fixed entities on right hand side of Eq. ( 133|) . 

In the absence of the electric quadrupole field, the s/, s„, Sj may be replaced 
respectively by appropriate M/,M„,Mj which are determined by the Kro- 
necker 5-function in Eq. ([7]), when the magnetic field B alone is present and 
is directed along {9b, (Pb)- Thus, cS and ctl, are replaced respectively by 
DLj. r i(pB,OB,0)* and (0b, ^b, 0) with Mj and Mi being fixed by left 

hand side. It may be noted that 0„ depends on M„ and the summation over 
n includes J2j„ Y.m„=-j„, with Q^" replaced now by Q^'^" whose elements are 
given by 

= dI,,mM, 0) z^^, M„(0, o)* . (35) 
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If the magnetic field B is along the Z-axis of ARF itself, the s/,s„,Sj in 
Eq. (!32l) may respectively be replaced by m/,m„,mj determined by the Kro- 
necker 5- function in Eq. Thus, cj, and c% are replaced respectively by 
^m'.mf and Sm'-rui, where rrif and rrii are fixed by left hand side of Eq. fl32|) . 
Therefore, the summation with respect to m'j and on right hand side of 
Eq. ( 132|) drops after making the replacements = 1 and 0% = 1. The ipn 

f i 

depends on m„ and the summation over n includes J2j„ J2m„=-Jn^ with Q^" 
replaced by whose elements are given by 

^m'lm'„ = ^m'^mn Sm'„m„ , (36) 

i.e., Q'^'^ gets replaced by a diagonal matrix with zeros everywhere except 
= 1 in Eq. il. 

It may be noted that the atomic transitions from ipi to ipn following absorption 
of uj' and from ipn^oipf consequent to the emission of uj are virtual transitions, 
which do not satisfy the celebrated Bohr criteria. This is in contrast to ab- 
sorption or emission represented by the = 1 term. They are real transitions 
which satisfy the Bohr criteria as already pointed out. The summation over n 
includes all atomic states ipn with different energy eigenvalues En- However, 
all of them do not contribute equally to Eq. (l32l) . The presence of 0„ on right 
hand side of Eq. (1321) indicates that one has to pay more attention to con- 
tributions coming from those states ipn with En close to Ei + uj' = Ef + uj. 
If there is an En such that En = E^ + u' = Ef + u, the contribution from 
this state alone overshadows all other contributions. The scattering is then 
referred to as resonance scattering. In particular, if Ei = Ef the terminology 
'two level resonance scattering' is employed. This is shown as (a) in Fig. [21 
where uj' = u. On the other hand if Ef > Ei as in (6) of Fig. O the resonance 
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scattering is referred to as three level resonance or fluorescence. 

If there is no electric quadrupole field and the atom is exposed only to a pure 
magnetic field B, such that the (2J„ + 1) states |J„M„) refer to distinctly sep- 
arated energy levels as in Zeeman effect, one can envisage resonance scattering 
taking place individually with each one of these taking the role of the upper 
level, as shown in (a) and (6) of Fig. [2|, if the condition Em^ = Ei + oj' = Ej + uj 
is satisfied. On the other hand, if gB < r„, the levels are not distinct and all 
of them contribute coherently to form a single line. This is referred to as quan- 
tum interference in the context of Hanle scattering, which is shown as (c) in 
Fig. [21 In contrast to Hanle effect where interference occurs between magnetic 
substates with the same J„, interference effects between states with different 
Jn have also been observed in polarization studies of Solar Ca II H-K and Na 
I Di and D2 lines [13], wherein it is mentioned that this can take place even 
when the lines are 3.5 nm apart. The general terminology, 'Raman scatter- 
ing' has been employed p^p3] to denote scattering, where contributions from 
several intermediate states are involved. In general, therefore, we may rewrite 
Eq. fl2T|) in the form 

r;,(k,/i;k',/i') = (^/|^(k,/x)|^„) , (37) 
where {ipv) represents a virtual state defined by 

1^.) = E<l^n); K) = 0„(^J^(k',/x')|^.) , (38) 

n 

which is clearly not an eigenstate of energy. In Raman effect, shown as (d) in 
Fig. [21 the lines corresponding to Ef < Ei are referred to as anti-Stokes lines, 
in contrast to those with Ef > Ei referred to as Stokes lines. 
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2-4 The dipole approximation 

If we neglect the spin-dependent second term in Eq. (1251) and employ dipole 
approximation e*^"" ^ 1 in Ak/x(r) given by Eq. (fT3l) . then we may express 
Eq. ([21D as 

{iju\A{k, ^ (^Je^ • p\iJi) ^ {iJi\S{k,fi)\ij^)* . (39) 

In the above equation the momentum operator p = — iv be replaced by 
[r, Hq] , to obtain 

{^P^\e^ ■ [r,i7o]|^z) ~ (^z - ^.)(^/'J^z), (40) 

if Eu, El denote the energy eigenvalues of ipu, i^i when considered as eigenstates 
of Eq. ([1]). We, thus, realize the Kramers-Heisenberg form represented by 
Eq. ^ of [H. 

3 The scattering matrix for atoms in external electric quadrupole 
and uniform magnetic fields 

The central result of the previous section is the derivation of the general 
expression for the on-energy-shell T-matrix element Tyj(/x, //'). If the incident 
radiation is in a pure state 

^. = Ec;'^;', (41) 

the amplitude for detecting the scattered radiation in a pure state 

^/ = Ec;^M, (42) 
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is given by 

Td£f,e,)=J2cf; d,,Tf,{fi,fi') , (43) 
mm' 

where IcJ^/p = 1 = I]^|c^P- On the other hand, it is more convenient 
to employ the density matrix formahsm [37|24ll25j to describe the states of 
polarization of the incident and scattered radiation, as it is more general and 
can handle mixed states of polarization as well. 



3. 1 The density matrix for polarized radiation 



The density matrix p for polarized radiation may be written as 



P=2 



1 3 



(44) 



in terms of the well-known [12] Stokes parameters (/ = Sq,Q = Si,U 



S2,V = S-i) and Pauli matrices a2 



a1 and the unit 



matrix o"o whose rows and columns are labeled by the left and right circular 
polarization states |/i = ±1) of radiation. Clearly, 



Sp = tr{a;p), P = 0,l,2,3 (45) 

where tr denotes the trace or spur. A column vector S with elements Sp,p = 
0,1,2,3 is referred to as the Stokes vector for polarization. If we consider 
Eq. ( 12T|) as a 2 X 2 matrix T with elements T^^/ = T/j(/i,/i'), the density 
matrix p of scattered radiation is given by 

P = T p' Tt , (46) 
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where p' denotes the density matrix of polarized radiation incident on the 
atom. Using Eq. (l45ll . we have 

Sp = ^ E tr{a; Ta} T^)S',, , (47) 

^ p'=0 

for the Stokes parameters of the scattered radiation, in terms of the matrix T, 
its hermitian conjugate T"!" and the Stokes parameters S'pi characterizing the 
radiation incident on the atom. 

3.2 The scattering matrix 

If the Stokes vector S' with elements (/' = S'^^Q' = S[,U' = S'^^V = S'^), 
characterizes the radiation incident on the atom, the Stokes vector S charac- 
terizing the scattered radiation may be expressed as 

S = 7^ S' , (48) 

where the 4x4 matrix TZ is referred to as the scattering matrix. Comparison 
of Eqs. ( l47l) and ( l48l) readily identifies the elements of 71 as 

^pp' = 2 ^ {crp)^L"^lT^,^L'{(^J)^l'^,"'{T^)^,"'^l" , (49) 

mm' 
m"m"' 

where we may use Eq. (1321) for T^^/ and note that (T''")^/"^// = T*,,^,,,, for which 
we may use the complex conjugate of Eq. fl32|) . We may thus write 

^MM' = E E (^rl'/m'^^m'^m^ (/i, /^') , 
m',m' 

(Tt)^,.,„ = t;„^„, = E E C;4-^-}'K' (/^"' /^'")* , (50) 

n' m'im" 
f » 

where 
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J * 

since Q^"' is hermitian. Using the above in Eq. (H9|) . we have 



m'im" 
f » 



(51) 



/ » 
m''m" 



f -i- 



^'-m'.'-m" 



^t(k',;x'")6^^"'^^(k,/i") 



St "f 

C „C I 

mm!, J J 



(52) 



Following Eq. flMl) . we may define hermitian matrices Q'^^ and ^^-^ through 
their elements 



'm'.m" 



^m'.^m" ' 



(53) 
(54) 



so that we may rewrite Eq. fl52l) as 



Tr 



^(k, /i)6;^M(k', /i')6^^'^^(k', /i'")e^^"'^^(k, /)6^^ 



(55) 



where Tr = J2m'j, denotes the Trace or Spur of the (2 + 1) x (2 + 1) matrix 
within the square brackets, which is defined through matrix multiphcation of 
the eight matrices, each of which is well-defined through Eqs. fl28l) . fl29l) . fl35|) . 
fl53l) . fl5^ for any specified atomic transition from an initial state ipi with 
energy Ei and total angular momentum Jj to a final state ipj with energy 
E'j and total angular momentum Jj, when the atom is exposed to a combined 
external electric quadrupole field and a uniform magnetic field. It may be noted 
that Si and Sf are fixed and the summation over n, n' includes summation over 
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3.3 The particular case of resonance scattering via electric dipole transitions 
between Ji = Jj = and Jn = ^ 



In this important particular case, which has often been investigated in the 
presence of pure magnetic fields, it is clear that Q^^ = Q'^^ = 1 in Eq. (155|) and 
L = L' = 1 in Eqs. (1251) and (12^ . so that we may write the trace appearing 
in Eq. fl55|) as 



Tr 



m"'m"" 



E ^(k',/i')^nk',/x'") 



^t(k,/i")^(k,/i) 



(56) 



We may use Eq. (fT5]l . in combination with Eq. (125]) . to write 



\J,{uJ')\'^^'^^"'Dl,A<t>\ 0\ 0)^1 , 0)* 



(57) 



and Eq. f[T^ . in combination with Eq. ([SI]), to write 



£t(k,/i")^:(k,/i) 



Using Eq. for f J and f Q^^' J , we may attach c^, c^,'„ along 

with (cTp/) n' to Eq. (!57|) . while we may attach c^„c^'„„ along with (o"^)/^"^* 
to Eq. dSSD, so that 
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= E E E /p'(A.,/..)F„„,(A„,ma)D^:^j0',e',O) , (59) 

Aa=0 Ha=—^a Tna=—^a 



Y: C.C„,,D^™^„(0,^,O)i^^„,(0,^,O)* 

= E E E /p(Ae,/ie)Fn,n'(Ae,me)*L>^V^(0,^,O). 

Ae=0 lXe=—\e me=—\e 



(60) 



Thus, we have 



1, 



2 2 



'^pp'^-^\Jl{^)?\Jl{^')?Yl,^ri<Pn' E E -^",n'(Aa,"^a)-^n,n'(Ae 

nn' Aa=OAe=0 

/,(Ae, /Xe)/p'(A„, /Xa)I)t^,(0, ^, 0)/?tM.('^'' 0) ' 



(61) 



where 



F„,„,(A,m) = EC(1,1, A;m;,-m:',m)(-l)<'c5,4"^, , (62) 
/,(A,/x) = ^C(l,l,A;/x', V",/^)(-1)"^'W"(^;)mV" ■ (63) 

A 

The /p(A, //) for p = 0, 1, 2, 3 may exphcitly be written as 



/o(A,/^) = ;^ 

/i(A,//) = -5a,2 



5a,0 + ^5a,2 



5m,2 - <^M,-2 



/3(A,//) = V^ 5a,i5^,0 ■ 



(64) 



It is interesting to note that F„^„/(A,m), in the particular case of an atom 
exposed to a pure magnetic field B directed along {9b, (ps) niay be written as 



Fn,n' (A, m) = Fm^^mL (A, 1^) 

= (_l)M;-rn ^(^^ ^. ^ ^ 

^AMa(6'b, 0_b) , 



(65) 
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in terms of the spherical harmonics. From Eq. (16T!) we can recover the Hanle- 
Zeeman scattering matrix of Stenflo [H] by setting the electric field strength to 
zero, and employing Eq. (1651) . The profile functions appearing in Eq. fl6ip are in 
the atomic frame. To obtain the Hanle-Zeeman scattering matrix, one needs to 
follow exactly the procedure outlined in ^14j , to transform 0„ to the laboratory 
frame through a Doppler convolution. In the weak field limit, we recover the 
well known Hanle scattering phase matrix of Landi Degl'Innocenti and Landi 
Degl'Innocenti [Ml, by assuming 4>m„<P*m' be independent of Mn,M^. In 
the strong field limit, gB is large compared to the line widths. Hence the 
three Zeeman component lines are well separated. If we neglect the coupling 
between the Zeeman substates (drop the summation over n in Eq. fl6T]) . by 
setting n = n'), we recover the restrictive phase matrix of Obridko [SH], which 
is basically a modified resonance scattering by individual Zeeman components. 

4 Numerical Results and Discussions 

The calculations presented in this paper are applicable to magnetic fields of 
arbitrary strength, and also the presence of quadrupole electric fields sur- 
rounding the radiating atom. To check the correctness of our derivation, we 
have reproduced the results of Stenflo (1998, Fig. 3) for the particular case of 
Hanle-Zeeman effect |40j|. In weak magnetic fields, pure Hanle effect prevails. 
In strong fields, the Zeeman effect is the dominant process. In intermediate 
fields, there is a smooth transition from weak field Hanle effect to the strong 
field Zeeman effect. These two effects exhibit relative dominance in differ- 
ent regimes of field strength, but they fundamentally overlap over the entire 
regime. 
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We consider the simplest case ofaJ = 0— s>0 type transition which 
produces a standard Zeeman triplet. In this Section we present the results of 
a single scattering experiment (see Eq. HSi) . We consider a 90° scattering of 
an unpolarized beam of radiation incident on the atom. The incident radia- 
tion is also assumed to be frequency independent (broadband pulse). Since 
S' = (1, 0, 0, O)""", the scattered Stokes intensity components are nothing but 
(/, Q, U, V) = (TZu, 7^21, 7^31, T^Ai), which measure the maximum degree 
of anisotropy for a given angle of scattering. The external magnetic field is 
assumed to be oriented along the Z— axis of the astrophysical (laboratory) 
reference frame (see Fig. [1]). 

The scattering is assumed to be frequency coherent in the laboratory frame. 
A Voigt profile function with a damping parameter a = 0.004 in units of 
the Doppler width Aud = {pqI c)\j2kT / Ma is employed to compute all the 
results. In Sect. 4.1 we present the results for pure magnetic case. In Sect. 4.2 
the results of scattering on atom immersed in pure electric quadrupole field 
is considered. Finally in Sect. 4.3 we consider the combined case of uniform 
magnetic and quadrupole electric fields. 

4-1 The pure magnetic field case 

Fig. [3] shows singly scattered polarization profiles for this case. The field 
strength B is chosen to represent the entire regime of Hanle-Zeeman effect 
{vb = 0.0008 - 2.5 in steps of a factor 5). The splitting parameter vb is de- 
fined as Vb = I'll where = eB/Airmc is the Larmor frequency. The 
geometry of scattering chosen by us (^' = 0°, 0' = 0°; 6 = 90°, (f) = 45°) 
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corresponds to Stokes / and Q given by 



/ = - Ml + ; Q = -i. (66) 

o '- -' 

The profile functions 0i and are now in tlie laboratory frame. They are 
obtained by a convolution of the atomic frame Lorentzian (see Eq. f l22]) ) with 
the Doppler profile. The real part of 0i for example, is a Voigt function whereas 
the imaginary part is a Faraday- Voigt function (see [H]). In the laboratory 
frame the frequency is expressed in dimensionless units (v = (z/q — v)/ ^vd)- 
Clearly, for weak fields the Zeeman splitting is not complete. Hence the / 
profile simply broadens without exhibiting a separation of the components. 
This kind of intensity profiles are very typical of Hanle and Hanle-Zeeman 
regime in the second solar spectrum of the Sun (see [12])- The case of vb = 2.5 
represents a strong field Zeeman effect, and we clearly see a well separated 
doublet. Since the line of sight (the scattered ray) is perpendicular to the 
magnetic field, according to conventional Zeeman effect theory (Zeeman effect 
treated as absorption/emission), one expects a triplet pattern in / profile, and 
a Q profile with vr component having opposite polarization compared to the 
two 0" components, along with U = V = Q. However, we now observe only 
a doublet in the / profile (see Fig. [3]), showing that the mechanism involved 
is indeed a 'scattering' process and not 'absorption followed by uncorrelated 
emission' process. With the help of classical theory of dipole scattering, one 
can argue that, for 90° scattering and for a magnetic field along the Z - 
axis, the vr component is not excited at all by the incident radiation. Only 
the components with electric vibration perpendicular to the scattering plane 
(containing the incident and scattered ray) are excited, and hence the two a 
components appear in the / profile. 
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From Eq. (166|) and also from the Fig. [3], we see that Q/I = —1, i.e., inde- 
pendent of frequency as well as the field strength. In other words scattered 
polarization is same as the well known non-magnetic pure Rayleigh scattering 
polarization. This is to be expected, because in the weak field regime, the Hanle 
effect is absent for vertical magnetic fields, and the scattered polarization can 
arise only due to Rayleigh scattering process. In the strong field regime, both 
the sigma components have linear polarization of equal magnitude, and scat- 
ter independently (see Eq. (166|) ). Therefore we obtain a maximum degree of 
linear polarization, namely Q/I = —1. We refer to this 'Zeeman scat- 

tering' or equivalently 'Rayleigh scattering in strong magnetic fields' (see also 
[TT1I23] ). 



4.2 The case of pure electric quadrupole field 

Fig. m shows the singly scattered Stokes profiles in the presence of a pure elec- 
tric quadrupole field surrounding the atom. The single scattering experiment 
is considered as in the pure magnetic field case. The electric field is charac- 
terized by two parameters, namely the electric splitting parameter A which is 
taken in the same way as magnetic splitting parameter, namely va = 0.0008 
- 2.5 in steps of 5, where va is electric splitting parameter in Doppler width 
units. The asymmetry parameter 77 is fixed at 1. The upper J„ = 1 level splits 
into 3 levels with energies — 2A, and +2A corresponding to the eigenstates 
V'l = |1, 0), ^2 = (|1, -1) - |1, l))/v^, and = (|1, -1) + |1, 1))/V2, respec- 
tively (see Fig. 2c of [21]). For the scattering geometry employed by us, the 
Stokes parameters are given by the analytic expressions 
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U = ^ 0203 + 030: 



* 

'2 ; 



V 




3.r 



(67) 



The profile functions 02 and 03 correspond to tlie eigenstates ijj2 and -03 re- 
spectively, and are given by 



in tlie laboratory frame. Here H and F are the well known Voigt and Faraday- 
Voigt functions [14j. From Eq. fl67|) . we note that U and V are generated purely 
due to the coupling between the eigenstates -02 and ip^. 

For weak electric fields {va < 0.1) the shapes of the / profile are not affected 
significantly, when compared to the corresponding pure magnetic case (see 
Fig. [3]). For va > 0.1 there is a blue shift compared to the pure magnetic 
case, which increases gradually for larger values of va until we get a doublet 
that is asymmetrically placed about the line center. The Q/I profiles have an 
interesting shape, in the sense that they are similar to the V/I profiles of the 
pure magnetic field case. For va = 0.0008 and 0.004, Q/I is extremely small. 
As Va increases, Q/I gradually increases, and the zero cross-over point of Q/I 
shifts toward the blue (see for eg. dash-triple dotted line in Fig. Hj). [/// and 
V/I show interesting behavior. Unlike Q/I, the ratio U/I takes largest value 
for Va = 0.0008, and is entirely positive (see upper most solid line in U/I panel 
of Fig. H]). For Va = 0.004 (dotted line) U/I takes both positive and negative 
values. As va increases we note that U/I goes to zero in the line core (except 
for = 2.5 - dash-triple dotted line which shows a small positive peak around 
V = 2.5), and becomes entirely negative and nearly constant in the wings. The 
shape of V/ 1 is similar to that of U/I (namely, taking constant values at line 



'2,3 — 



H{v V)VA, a) + 2 z F(v + (1 T r])vA, a), 



(68) 
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core and line wings with a swift transition around v ~ 3). The magnitude of 
V/I initially increases with va and then rapidly decreases toward zero, as va 
increases. For va > 0.1, the V/I is nearly zero throughout the line profile (see 

Fig.SD. 

If the mechanism involved is a pure emission process, then one would expect 
all the three wave functions to contribute, and produce line components 
in Stokes /, along with Q = \/ = 0, and f/ ^ (see Eqs. (69)-(72) in p¥]). 
However when the interaction of radiation is treated as scattering (represented 
through angular correlations between incident and scattered ray), we see that 
only a doublet is seen in Stokes / (since ipi is not excited at all according to 
Eq. (!67|) ) with Q, U and V non-zero. This is the essential difference between 
the spontaneous emission process and the scattering process in quadrupole 
electric fields. We have also computed the Stokes profiles for the rj = 0.5 case, 
and find that they do not differ qualitatively from 77 = 1 case, except for 
changes caused by different amount of level splitting. 

4-3 The case of combined magnetic and quadrupole electric fields 

In Fig. [5] we show the Stokes profiles for this case. We employ ratio A/B = 
0.5 (which defines the relative strength of electric field with respect to the 
magnetic field), and the asymmetry parameter rj = 1. The splitting parameter 
vb (= 0.0008 - 2.5) is employed as in the pure magnetic case. The simplest 
geometry of the combined magnetic and quadrupole fields with B along the 
Z-axis of the PAF, and PAF itself coinciding with the ARF, is employed in the 
computation of the results in this section. The scattering geometry is same as 
in Figs. [3]and|H The upper level = 1 is split into three levels n = 1,2,3 
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with energies 



El = -2rB ■ ^2,3 = [rT (rV + 1)'^']^, (69) 



where r = A/B. The corresponding eigenstates are given by (see p 



iJi = |1, 0) ; ij2 = h\l, -1) + 62II, 1) ; ^3 = -&2I1, -1) + h\l, 1), (70) 



where the interference coefficients are defined by 



rr] 



+ l + (rV + l)'/' 



-1 1/2 • 



rT]-l + {r'^rf + 1)^/^ 

y,2yyl _j_ ]^ _|_ y,rq^rp2rq2 _j_ ][^l/2 



1/2- 



(71) 



Fig. [6] shows the energy level splitting for spin-1 upper level, exposed to the si- 
multaneous presence of an external electric quadrupole and uniform magnetic 
fields. Notice that "level-crossing" occurs when the electric field strength in- 
creases. It is worth noting that these levels are not pure states, but superposed 
states. The cross-over occurs for only a single value of electric quadrupole field 
strength at which the levels become degenerate. 



Fig. [5] presented in this paper can be interpreted using the panel (c) of Fig. [61 
For r = 0.5, the energy eigenstates Ei and E2 are below the A = B = Q 
reference line. For the geometry chosen, the state ipi = \1, 0) is not excited, as 
can be seen from the following analytic expressions for the scattered Stokes 
parameters : 
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(72) 



The profile functions 02 and 03 in tlie above expressions correspond to tlie 
eigenstates -02 and -03 respectively, and are given by 

02,3 = H(Y+[rT^r^V^ + l)vB,a^+2tF(Y+{rT^JrW + '^)vB,ay (73) 

in tlie laboratory frame for any rj. 

The Stokes I profiles in the combined case are quite similar in amplitude 
and shape to those in the pure magnetic case, except for the position of the 
component lines (for vb > 0.1). For vb = 0.0008 and 0.004 (solid and dotted 
lines which nearly overlap on each other m Q/I panel of Fig. [5]), the Q/I 
profiles are very different compared to the corresponding pure electric field 
case, and are entirely negative. For vb > 0.004, these Q/I profiles resemble 
their corresponding counterparts of the pure electric field case in terms of the 
frequency dependence. Clearly, there is a decrease in the magnitude of Q/I 
when compared to both the pure magnetic and the pure electric case. 

U / 1 m the combined case is quite similar to the pure electric case (compare 
U / 1 in Eq. fl67|) with Eq. fl72|) for r = 0.5 and 7] = 1). In other words, we can 
conclude that for this particular geometry the entire frequency dependence of 
the U/I comes from the electric field effect. The only difference is in vb = 2.5 
(dash-triple dotted line) case, where in the small peak is enhanced as well as 
it is now centered around 1.5, unlike in pure electric field case. V/I profiles 
in combined case bear resemblance to Q/I profiles in shape, differing only in 
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magnitude and sign (for vb > 0.0004). Clearly, the quadrupolar electric field 
produces an additional depolarization in Q/I (compared to the Q/I of pure 
magnetic case in Fig. [3]), and a rotation of the plane of polarization (or gen- 
eration of new U/I). Therefore in the special geometry of vertical magnetic 
fields (B parallel to the Z-axis of PAF and ARE simultaneously), if one ob- 
serves strong U/I and V/I signals as well as relatively smaller Q/I signals in 
a 90° scattering (eg. as in the case of extreme limb observations of the Solar 
chromosphere; or like the special scattering geometry that we considered for 
discussion), it could indicate the asymmetries arising from the quadrupolar 
electric fields surrounding the atom. Notice that in the absence of electric 
fields, U/I = V/I = and Q/I = —1 for the vertical magnetic fields (see 
Fig. [3]). Another diagnostic indicator of quadrupolar electric fields is the net 
blue shift of the Stokes profiles, unlike the linear Stark effect which produces 
symmetric shifts with respect to the line center. 

For the sake of discussion, we have computed the scattered Stokes profiles 
for rj = also in the combined case. This case is interesting, because, from 
Fig. E] we observe that for r < 0.5, the splitting is independent of rj, namely 
the splitting pattern is same for both rj = and 1. As a result, the Stokes / 
for both 1] = and 1 cases are nearly identical. This can be understood from 
Eq. ( 172]) by setting r] = 0. The profile functions 02,3 are now given by Eq. (1751) 
with rj = 0. We note that the real part of 02 for t] = and 1 respectively 
are iJ(v — 0.5vb, cl) and if(v — 0.62vb, cl)- Similarly for 03, corresponding real 
parts are II{v + 1.5vb, cl) and II(y + 1.62vb, cl). For this reason, the Stokes / 
as well as U/I profile for 77 = and 1 are nearly identical. However the Q/I 
and V/I profiles are quite different (see Eqs. (!72|) ). The Q/I profiles of rj = 
case have shapes and magnitudes similar to the V/ 1 profiles of Fig. |H The 



30 



V 1 1 profiles of 77 = case liave sliapes and magnitudes similar to Q/I profiles 
of Fig. m except for a sign difference. 

5 Conclusions 

The scattering matrices for the combined effect of electric quadrupole field and 
uniform magnetic fields (of arbitrary strength) are derived using quantum elec- 
trodynamic approach. The scattering matrix for Hanle-Zeeman effect is vali- 
dated by comparing with the published results of Stenflo [14j . The quadrupole 
electric field is characterized by strength and asymmetry parameters, which 
produce unique diagnostic signatures that may be employed to detect the elec- 
tric charge distribution asymmetries in the solar atmosphere. The theoretical 
formulation is quite general and can handle not only the simplest case of a 
triplet ( J = ^ 1 ^ 0) transition that is employed for illustrations in this pa- 
per, but also an arbitrary choice of quantum numbers. We have demonstrated 
the properties of the coherence or interference phenomena, like strong field 
Zeeman scattering, Hanle and magnetic Raman scattering that are important 
in the interpretation of spectral lines in the second solar spectrum. 
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Fig. 1. The scattering geometry: (Xq, Yg, Zq) refers to the Principal Axes Frame 
(PAF) characterizing the electric quadrupole field. The radiation is incident along 
{9' ,(f)') and scattered along {6,^) with respect to the astrophysical reference frame 
(ARF) denoted by (X,Y, Z). The magnetic field B is oriented along {9b-,^b) with 
reference to PAF and {9b,4'b) with reference to the astrophysical reference frame 
(the azimuthal angles and (f)B are not marked in the figure). 
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Fig. 2. Level diagrams sliowing the atom-radiation interaction processes discussed 
in this paper, (a) two-level resonance scattering process, (b) three-level fluorescence 
scattering process, (c) Hanle scattering process in weak magnetic fields and (d) the 
general case of Raman scattering. 
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Fig. 3. Stokes Profiles for Pure magnetic field case. Magnetic field is oriented 
along the Z -axis of the ARF (see Fig. [1]). The model parameter employed are: 
a = 0.004, the scattering geometry defined by 6*' = 0°, 0' = 0°; 61 = 90°, 4) = 45°, 
and the magnetic field strength defined by vb = 0.0008 (solid line), 0.004 (dot- 
ted line), 0.02 (dashed line), 0.1 (dash-dotted line), 0.5 (dash-double dotted line), 
and, 2.5 (dash-triple dotted line). Notice a constant degree of linear polarization 
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100%. 
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Fig. 4. Stokes line profiles for pure quadrupolar electric field case. The 
model parameters employed are: a = 0.004, with the scattering geometry 
61' = 0°, ^' = 0°; 6* = 90°, = 45°, and the electric field strength defined through 
the splitting parameter va = 0.0008 (solid line), 0.004 (dotted line), 0.02 (dashed 
line), 0.1 (dash-dotted line), 0.5 (dash-double dotted line), and, 2.5 (dash-triple dot- 
ted line), and the asymmetry parameter jy = 1. Notice the blue shift of the profiles 
compared to the pure magnetic field case for large values of va (> 0.1). 
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Fig. 5. Stokes Profiles for combined quadrupolar electric and uniform magnetic field 
case. The ratio r = A/B = 0.5 and the asymmetry parameter r/ = 1. Other input 
model parameters, scattering geometry and the line types are same as in Fig. [3l 
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Fig. 6. Energy level diagrams showing the combined effect of magnetic and electric 
quadrupole fields. The panels (a, b, c) represent different asymmetry parameters r\. 
The energies Ei, E2, and £^3 corresponding to the states and ^3 respectively 

are plotted. The value r = corresponds to the pure Zeeman case. 
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